A TFD model for the Electrospheres of Bare Strange Quark Stars 
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We study the layer of electrons on bare strange star surfaces, taking the Dirac exchange- 
energy into account. Because electrons are fermions, the electron wave function must be of 
exchange-antisymmetry. The Dirac exchange-energy originates, consequently, from the exchange- 
antisymmetry of electron wave functions. This consideration may result in changing the electron 
distribution and the electric field on the surface of bare strange star. The strong magnetic field 
effect on the structures of the electrospheres is also discussed. 
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We are unfortunately not certain about the real state 
of matter in pulsar-like stars even 40 years after the dis- 
covery of pulsar. Nuclear matter (neutron stars) is one of 
the speculations.while quark matter (quark stars) is an 
alternative [e.g. , B, B 0, II B @, 0, B ©, \m ■ The observa- 
tional features could be very similar for normal neutron 
stars and crusted strange quark stars, but the bare quark 
surface should be very useful for us to distinguish those 
two kinds of compact objects It is therefore worth 
studying the structure of electrospheres of bare strange 
quark stars in details [for the previous researches, see, 
e.g.M- 

If the Witten's conjecture [l3[ (the strange quark mat- 
ter in bulk is the most stable strong interaction system) 
is correct, strange quark star with almost equal num- 
bers of it, d, and s quarks would exist in case that the 
baryon density is higher than the nuclear density. In the 
strange quark matter, the three flavors of quarks may 
keep a chemical equilibrium, and the whole of the quarks 
will be electropositive because of massive s quarks. In 
order to make sure the strange quark matter is globally 
electroneutral, there must exist proper quantity of elec- 
trons. 

Alcock et al. [|| thought that (1) a strange quark star 
would accrete inter-stellar matter, and then a crust would 
form, wrapping the strange quark star, and (2) a bare 
strange star can not manifest itself as a radio pulsar 
because of being unable to generate a magnetosphere. 
However, it is argued [H, EH, that there are some ad- 
vantages in understanding the observational features if 
radio pulsars are bare strange stars and that the circum- 
stellar matter could not be accreted onto the star's sur- 
face at all (i.e., the quark surface keeps bare). The elec- 
trons can only extend thousands of femto-meters above 
the quark surface, and an extremely strong electric field 
E~5x 10 17 V/cm forms @, [H, QJ, El, beinghigher 
than the critical field E CT ~ 1.3 x 10 16 V/cm [20[. It 
is worth noting that the Usov mechanism (e ± -pair pro- 
duction in degenerate electron gas with strong electric 
field) [21], |22j, |23j may dominate in the thermal emission 
of a bare strange star with a surface temperature T s , 
6 x 10 8 K < T s < 5 x 10 10 K. 



Sanwal et al. [24] discovered two absorption features in 
IE 1207-5209's X-ray spectrum^OJ keV, 1.4 keV. Based 
on the works of Ruder et al. [25| and Mori and Hailey 
et al. [26(, Sanwal et al. thought the two absorption 
features originated from the atomic transitions of once- 
ionized helium in the neutron star atmosphere with a 
strong magnetic field. They do not think it is likely to 
interpret the two absorption features in terms of electron 
cyclotron lines for four points of view. Howsoever, these 
points were criticized by Xu et al. [27[, who suggested 
that the two absorption features could be of electron cy- 
clotron lines from the bare strange quark star surface. 
Xu et al. proposed a model with a debris disc around a 
bare strange quark star, which could slow down the spin 
and decrease the inferred strength of the dipole magnetic 
field so as to fit the energy differences of Landau levels 
(~ 0.7 keV). It is then very necessary for us to under- 
stand the exact electron distribution on the bare strange 
star surface in terms of cyclotron lines. 

The extreme relativistic electron distribution of the 
bare strange star surface has been calculated via Thomas- 
Fermi Model [1, E3 , but the Dirac exchange-energy [28| 
was not taken into account up-to-now. The wave func- 
tion must be of exchange-antisymmetry for fermions, and 
the Dirac exchange-energy originates consequently. For 
the Dirac exchange-energy plays an important role in the 
low-density electrons and make a dramatic effect on the 
electron distribution, we need to take it into account, 
applying the Thomas-Fermi-Dirac (TFD) model for bare 
strange stars. 

The theoretical basis of the Thomas-Fermi-Dirac 
Model is the density functional theory, called DFT for 
short. Via this theory, we could obtain the Thomas- 
Fermi-Dirac Equation of the electron density. 

In DFT, the ground state energy of the system is the 
functional of the ground state electron density, which 

is m 



E Q = E [n (r)] 



(1) 



where Eq is the ground state energy of the system, hq 
is the ground state electron density. The ground state 
energy of the electron gas is composed of three parts (we 
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used the different form from ref. 



Eo[n {r)} = -e \ n (r)V {f f ) df 

s 2 

\n-r 3l 



(2) 



where |\&o)is the ground state vector of the electron 
gas, the first term of right side is the outer poten- 
tial energy of the system, the second term is the inter- 
electronic Coulomb energy and the last one is the elec- 
trons kinetic energy. Here we can divide the inter- 
electronic Coulomb energy into two terms, the classical 
inter-electronic Coulomb energy 



no(f)n (f') 



drdf' 



and the inter-electronic exchange-energy, also called 
Dirac exchange- energy, E ex . And then we can put the 
terms of the outer potential energy and the classical inter- 
electronic Coulomb energy together 

— e J no(r)V(r)dr : 

while putting the terms of the Dirac exchange-energy and 
the electron kinetic energy together 



G[n (r)]dr, 



where G[no(f)] is also the functional of the ground state 
electron density no(r). So we get the expression 

Eo[no(r)] = -e I n (r)V(r)df 



G[n (r)]dr. 



(3) 



For the ground state energy is the minimum of the sys- 
tem, it must be content with the variational equation [28] 

SE [n (f)} = -e J 6n (r)V(r)dr + J SG[n (r)]df 

= -e [ 5n (r)V(r)df + [ SG P l °^ Sn (f)df= 0. (4) 
J J 5n (r) 

As the electron number is conservative, it must be con- 
tent with the variation equation [28| 



SN = / 5n (r)df = 0, 



where N is the total number of electron. If we take a 
Lagrange multiplier [28| . — p, called chemical potential, 
we can joint the two variational 0-equations above into 
one variational differential equation, 



5n {r) 



(6) 



In classical electrodynamics, 
reads, 

\7 2 V = -47rp e 



the Poisson equation 



(7) 



where V is the potential produced by the electrons and 
the quarks, while p e is the charge density of the elec- 
trons and the quarks. However, out of the strange quark 
matter surface, there are no quarks. If we take the vari- 
ational differential equation into the Poisson equation, 
we'll obtain the Thomas-Fermi-Dirac equation, 



V 



2 8G[no(r)} 
Sn (r) 



47re 2 (n (r) - n q ), 



(8) 



where e(n (r) — n q ) = —p e , en q is the charge density of 
quark. For simplicity, we would assume that the strange 
quark matter is uniform global electropositive, and n q is 
then a constant n„ = n„„ in the strange quark matter, 



~ lb qp 

while out of the strange quark matter surface, n q = 0. 
For the thickness of the electron gas out of the quark 
matter surface is so thin to the radius of the strange 
quark star that we can ignore the curvature of the star 
surface, we can take ^ instead of V. If we are able to 
obtain the apparent relation between G[no(f)] and no(r), 
we can take d instead of 5. So we get the new-formed 
TFD equation 



d \ dG(n ) 
dz J dno 



4ire 2 {n - n q ), 



(9) 



where z is the height to the quark matter surface. When 
z > 0, it is out of the surface, z < 0, it in the strange 
quark matter, z = 0, it on the surface. If the both sides 

of the equation above are multiplied by d( d ^"°-* )(see 
ref. [28j]), it can be transformed into a lst-order ordinary 
differential equation 



d dG(no) 



d ( 4?re ( (n - n q ) 



dG(n ) 



G(n ; 



N dz dno J J V V driQ 

(10) 

It can be inferred that at the infinite height, n — and 
dno/dz = 0, so out of the strange quark star surface, the 
equation above can be integrated from the infinite height, 
simplified to the expression (for z > 0) 



d dG(rto) y 
dz dno J 



Aire I no 



dG(n ) 
dn 



G(no) . (11) 



In the equation above, the left side of the equal sign is 
a square, in order to make sure the equation has a real 
number solution, the right side must be lager than or 



(5) 

equal to zero, i.e. 



n - 



d G(n ) 
dn 



- G{n ) ^ 0. 



(12) 



If we obtain the exact relation between G(no) and no, 
we can calculate the minimum of no via the in-equation 
above. With the increase of the height, no decreases con- 
tinuously to the point of its minimum and, after reaching 
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the point, will jump to zero directly for appropriating to 
the boundary condition. 

Deep in the strange quark matter, z — > — oo, it can 
be inferred that the electron charge density equals to the 
three flavors of quarks charge density [l4[ , 



no(z = — oo) = n q . 
If we integrate the equ.(17) from z = 

dG(no) 



(13) 



-oo,we '11 get 



1 / d dG(rao) 



: 47re z 



, , .„ , ■ (G(n ) - G(n qp )) 

az dno / V dno 

(14) 

On the strange quark matter surface, z = 0, the elec- 
tron density and the gradient of the electron density are 
continuous, so from the equ.(ll) and equ.(14), we get the 
equation below 



dG(n ) 



dno 



■ G(n ) = (no - n qv ) 



dG(no) 



(G(n ) 



dG(n ) 
dn 



+ G(n qp ), 



dG(n ) I 
dn 



G(n q 



■ G(n 9p ))(15) 
(16) 
(17) 



Consequently, if we obtain the expression of G(no) and 
the value of n qp , we could work out tiq via the equ.(17) 
as the boundary condition of the equ.(ll). Howsoever, 
we describe the strange quark matter by its chemical po- 
tential that is also called Fermi energy, 



Hc(3Tr 2 n qp ) 3 , 



(18) 



so we could work out n qp from fi q . 

Usually, the local density approach (LDA) is an effec- 
tive way to describe the Dirac exchange-energy. In this 
way, the Dirac exchange-energy reads [291 ] 



-Bex (n ) 



-2e 2 



3n 
8tt 



(19) 



Note that G(n ) is composed of the Dirac exchange- 
energy E ex (no) and the electron kinetic energy T(no). In 
the extreme relativistic case for electron, we could obtain 
the electrons kinetic energy via Fermi-Dirac statistics 



(T) = -hc(3ir 2 n )3 . 



(20) 



Actually, we've ignored the influence caused by tem- 
perature. While for the generally relativistic case, the 
means of the electrons kinetic energy could also be 
worked out via Fermi-Dirac statistics 



(T) = 



3 3 , 2 2 
+ -(37r 2 no) 3 

4 8 \m e c 



1 + (3ir 2 no)3 



(37r 2 n ) 1 (— — ) sinh 1 ^(37r 2 n )i (— " ~ ljm e c 2 (21) 



From the definition of G(uq), we get the exact relation 
among G(n ), E xc (n ) and (T)(n ) 



G(no) = n {(E xc (n Q ) + (T)(n )}. 



(22) 



For the extreme relativistic case, with the expression of 
G(n ), TFD differential equation and the typical chem- 
ical potential (e.g., 20 MeV [3]) of the quark in the 
strange quark matter, we could calculate the electron 
distribution out of the strange quark matter surface. If 
we ignore the Dirac exchange-energy, TFD model will de- 
generate to TF model calculating the extreme relativistic 
electron distribution out of the bare strange star surface. 
For the generally relativistic case, with the similar con- 
ditions, we could work out the distribution. The distri- 
butions of the two models TF and TFD of the two quark 
chemical potentials, 2 MeV and 20 MeV, are illustrated 
in the Fig.l. Also, we can work out the electric field, illus- 
trated in the Fig. 2. The maximum of it is ~ 10 17 V/cm. 
If we enlarge the range of the typical chemical poten- 
tial of the quark in the strange quark matter, choosing 
2 MeV, 20 MeV and 200 MeV separately, we could get 
three kinds of the distribution, illustrated in the Fig. 3. 

As in Fig.l, in the high-density domain, TFD- 
generally-relativistic electron density is about 2 orders of 
magnitude more than TF-extremely-relativistic electron 
density, while TFD-electron density is much less than 
TF-electron density in the low-density domain. The dif- 
ferences reflect the influence of Dirac exchange-energy. 
It is the Dirac exchange-energy that makes the electrons 
high above much lower, and it seems to put an attractive 
force. In Fig. 3, we observe that, in the low-density do- 
main, the electron distributions of different quark chem- 
ical potentials (i.e., V q ) have few differences. 

If the strong magnetic field is taken into account, the 
sum of the electron kinetic energy and magnetic energy 
is quantized (30| . 



E = yj m 2 c 4 + c 2 p 2 + rri^tvuj 



Bn, 



(23) 



here h = til + & + h, n L is the orbit quantum number, 
a = ±1/2 is the 3rd component quantum number of the 
electron spin, and cob = eB/mc is the cyclotron angular 
frequency in the magnetic field. With the increase of the 
height, the electron density n(z) will decrease to a critical 
density of n* [30| . 



2m 2 ix)BC [2huj 



K 2 



(24) 



The electrons above the height of n* will play an impor- 
tant part in the abstraction of X-ray [30( . If we take the 
standard strength 10 12 G of the dipole magnetic field, we 
can work out n* ~ 10~ 10 fm~ 3 . In Fig. 3, we can find 
out the point of n*/N = 10~ 2 . Only the electrons above 
the height of n* can absorb the X-rays emitted from the 
strange quark matter. Because TFD-electron density de- 
creases more rapidly after the n* than TF-electron den- 
sity, there are much fewer electrons which can absorb the 
X-rays. Therefore, one has to take the Dirac exchange- 
energy into account when one tries to study the electron 
distribution and the X-ray absorbtion features. 

In Fig.3, the point of n*/N = 10~ 2 for B - 10 12 G 
could be found and the three kinds of the electron dis- 
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tribution in low-density domain are similar, the distri- 
bution of the electrons, which are able to absorb the X- 
rays, is then not sensitive to the chemical potential of the 
strange quark matter (V q , for instance). Nevertheless, if 
the n* were larger, the distribution would become sensi- 
tive to the chemical potential of the strange matter. It 
is then possible that we might obtain the strange quark 
matter chemical potential by studying the X-ray absorb- 

3 3 

tion features. In equ.(24), n* — ^1^ 2 2 oc B? , we 

may need a strong dipole magnetic field, B > 10 13 G 
(n* /N = 10 _1 ). In one word, if we observe the bare 
strange star with the stronger magnetic field, the ab- 
sorption features of its X-ray spectrum may provide us 
a possibility to know its inner strange matter chemical 
potential. 

Summarily, the issues studied in the paper make it nec- 



essary for us to take the Dirac exchange-energy into ac- 
count when we consider the electrosphere of bare strange 
stars. The differences between the electron densities in 
TF and TFD models are significant. If the exact distri- 
bution of the electrons on the bare strange star surface 
is obtained and if the magnetic field of the bare strange 
star is strong (B > 10 13 G), one may probably obtain the 
electron chemical potential of the strange quark matter 
by investigating the thermal absorption spectra. 
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FIG. 1: This figure illustrates the differences of the electron 
density among the different models with different quark chem- 
ical potentials. N = {^f = 1.7 x 10~ 8 fm~ 3 . Z = ^ = 
3.9 x 10 2 fm. z is the height to the strange quark matter sur- 
face, n is the electron number density of z. In the extremely 
relativistic case, the difference between the electron density 
curves of TF and TFD models is not significant when the den- 
sity n is very high, but is significant when n is lower. In the 
high-density domain, TFD-generally-relativistic electron den- 
sity is about 2 orders of magnitude more than TF-extremely- 
relativistic electron density, while TFD-electron density is 
much less than TF-electron density in the low-density do- 
main. 
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FIG. 2: This figure illustrates the differences between the 
electric fields, E, in the TFD model and the TF model. The 
maximums of the electric field in the two models are both 
stronger than the critical field ~ I0 16 V/cm. Z = — — = 
3.9 x I0 2 fm. 




FIG. 3: This figure illustrates the differences of the electron 
density among the different quark chemical potential in the 
TFD model, z is the height to the strange quark matter sur- 
face, n is the electron number density of z. In the high-density 
domain, the electron distributions of different quark chemical 
potentials are very different. Nevertheless, in the low-density 
domain, the three electron distributions have few differences. 
The lowest horizontal line is correspondent to n*(B = 10 12 G). 
The middle horizontal line, which is to n*{B = 10 13 G), could 



tell the difference between V„ 



20 MeV and V a 



2 MeV. 



The highest horizontal line, which is to n*(B — 10 15 G), may 
tell the difference between V = 200 MeV and V„ = 20 MeV. 



